Manipulation of extreme events on scale-free networks 
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Extreme events taking place on networks are not uncommon. We show that it is possible to 
manipulate the extreme events occurrence probabilities and its distribution over the nodes on scale- 
free networks by tuning the nodal capacity. This can be used to reduce the number of extreme events 
occurrences on a network. However, monotonic nodal capacity enhancements, beyond a point, does 
not lead to any substantial reduction in the number of extreme events. We point out the practical 
implication of this result for network design in the context of reducing extreme events occurrences. 
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I. INTRODUCTION 

The study of dynamical processes [l| and extreme 
events 0, y| on complex networks has become an im- 
portant topic of research interest both for its inherent 
scientific understanding and possible applications. The 
main motivation being that many extreme events such 
as the floods, congestion in internet and other communi- 
cation networks, power black-outs and traffic jams take 
place on complex networks. In particular, the effect of 
these extreme events is tremendously high in terms of 
its impact on human lives, property and productivity 
This is evident, for instance, in the life and property lost 
due to floods and the working hours lost due to power 
black-outs and traffic jams [5|. It would be beneficial if 
some form of control or at least a possibility of manipu- 
lating the occurrence of extreme events can be achieved. 
In this work, we show that we can manipulate the ex- 
treme event occurrence probabilities on networks and we 
examine the extent to which this is feasible. 

We assume that a suitable model for transport is de- 
fined on any network such that the events in any of its 
node is proportional to the flux x{t) passing through it. 
For instance, flux could be the volume of water flowing 
in a river or the number of information packets passing 
through a router in an IP network. Extreme events are 
associated with exceedences of the flux above a certain 
threshold, i.e., an event is called an extreme event at time 
T \i x{t ~ T) > Xee where Xee is the threshold used to 
identify extreme events. In this setting, we focus on the 
extremes, arising due to inherent fluctuations, in the flux. 
Clearly, then the extreme events cannot be avoided alto- 
gether in any flnite system. They can at best be partially 
mitigated by tuning a suitable network parameter. This 
is the theme we explore in this work. As this approach 
allows some freedom to manipulate EE on networks, it 
is similar in spirit to the emerging interest in controlling 



the dynamics on complex networks [6|. 

Our transport model is the dynamics of random walk- 
ers on networks. Though random walk on networks were 
studied earlier [Tj, |8| , they were not focussed on the ex- 
treme event (EE) properties. In a recent work 0,13], the 
threshold for extreme events Xee was taken to be pro- 
portional to the typical size of the flux passing through 
a node, i.e, Xee ^ '^f — {.fV^^i where (/) and aj are 
the mean and the standard deviation of the flux passing 
through the node. In the present work, we choose the 
threshold to be 



cx a/ = (/)", a>0. 



(1) 
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Main motivation for this choice arises from the results in 
Ref. which show that cr / oc (/)" and in particular 
the value of a depends on factors such as the resolution 
of the flux measurements and the noise in the number of 
walkers. 

Another motivation arises because Xee can be inter- 
preted as the nodal capacity [ll| and extreme events oc- 
cur when flux exceeds the nodal capacity. The empiri- 
cal and modeling studies [l2[ on the relationship between 
mean load and mean capacity in a network do not address 
the question of extremes in load fluctuations. Extreme 
events result from large fluctuations in load and these 
are mostly responsible for temporary local network fail- 
ures. Our experience suggests that most of such extreme 
events, say the vehicles piling up at a traffic intersection, 
arise due to limited throughput capacity of the nodes. 
Then, one possible solution to reducing the extreme event 
occurrences would be to increase the handling capacity 
of the node. Guided by this intuition, we vary a which 
is a proxy for the nodal capacity such that larger values 
of a correspond to larger handling capacity of the node. 
In this work, we show that it is possible to manipulate 
the probabilities for the occurrence of extreme events on 
the nodes of a scale-free network by tuning a. 

A more significant result is that tuning a beyond a 
certain point does not lead to any significant reduction 
in the number of extreme events. This result has im- 
portant implications when we opt for capacity building 
as the route to mitigate extreme events on networks. 
Capacity addition invariably comes at a heavy price. 



2 



For instance, building bridges at a road intersection is 
one method of capacity addition that might increase the 
throughput across the node. In communication networks, 
additional servers and switches might be needed at a 
node to smoothly handle excess traffic flowing through 
it. All these interventions require expensive infrastruc- 
tural changes. In view of such high costs involved in 
this effort, it is important to ask if increase in capacity 
will lead to proportionate decrease in the likelihood of 
extreme events. The results in this work furthers our in- 
sight into this question in the context of extreme events 
on complex networks. 
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II. RANDOM WALKS AND EXTREME 
EVENTS ON NETWORKS 

We model the flux as a random walk process exe- 
cuted by W independent walkers on a fully-connected, 
unweighted network. In this case, the distribution of 
walkers passing through z-th node is given by 0, 0| , 

/.H = Q pt ii-Pir^'^- (2) 

The stationary probability pi to flnd a walker on i-th 
node is 



This measures the extent to which the walkers are at- 
tracted to i-th node and it depends on its degree. The 
flux of walkers being binomially distributed, the mean 
{fi) and variance af of the flux passing through i-th node 
can be immediately written down and both depend only 
on degree k^. 

In the spirit of extreme events statistics [9j , we define 
an event on node i to be extreme ii w > qi, the threshold 
Qi to be determined below. Then, the probability for the 
occurrence of extreme event on node i is given by 

w 

F^{q^)^ ^ Mw)^Ip{[q,\+l,W-[q,\), (4) 

where [.J is the floor function and Iz{a,b) the incom- 
plete Beta function 13]. We note that the extreme event 
probability depends on the choice of threshold qi and 
this serves as a handle to tune and control extreme event 
probability J-i. In this work, the threshold is chosen as 

g, = (/,)+wr«I^-^-fm(vF— ^) ,(5) 

where m,a > are real numbers. The magnitude of ex- 
treme event m scales the probability for extreme events 
and does not qualitatively change it. Hence, a pro- 
vides a handle to manipulate the extreme events by tun- 
ing the threshold qi. 



FIG. 1: The probability for the occurrence of extreme events 
as a function of degree k, shown as log-log plot, for a = 
0.40, 0.43, 0.47, 0.50, 0.52. Both the analytical (Eq. [4| and 
simulation results are shown. The shaded region on the left 
(right) corresponds to small degree nodes (hubs). 



A. Manipulation of extreme events 

To illustrate the idea of manipulating the extreme 
events, we display the probability for the occurrence of 
extreme events in Fig. [TJ The simulation results shown 
in this figure (and in the rest of the paper) are obtained 
from random walk simulations on a scale-free network 
(degree exponent : 7 = 2.2) of size N — 5000 and 
W = 39830 averaged over 100 realizations for each value 
of OL. Clearly, there is a good agreement between the an- 
alytical (Eq. HJ and simulation results. Note that as a 
is increased from 0.4 to 0.6, the probability for the oc- 
currence of extreme events Fi{k) changes systematically. 
For the hubs (fc > 80), Fi changes by nearly 2-3 orders 
in magnitude. A significant feature is that for a < 0.47, 
extreme event probability Fi is higher, on an average, 
for the hubs (fc > 80) when compared to the small de- 
gree nodes {k < 20). This feature reverses for a > 0.47, 
i.e., Fi is higher for small degree nodes compared to that 
of hubs. For a ~ 0.47, the probability Fi is approxi- 
mately independent of the degree of the node (ignoring 
the local fluctuations). By tuning a we obtain a range of 
behaviour for the EE occurrence probabilities. 

Physically, results in Fig [T] can be interpreted as fol- 
lows. Note that a = 1/2 and m = 1 in Eq. [5] represents 
the typical size Wtyp of flux through a node in the random 
walk environment. We can take the threshold qi to rep- 
resent the capacity of the node to handle this flux Wtyp. 
Then, a < 1/2 would imply Wtyp > qi, i.e, flux is more 
than the capacity of the node to handle it. This can be 
thought of as the congestion-like situation. On the other 
hand, a > 1/2 corresponds to wtyp < g^, in which typical 
size of flux is smaller than the capacity of the node and 
hence congestion-like situation is less likely to happen. 
Fig. [1] implies that as a is varied, the probability of EE 
on the hubs shows larger variability than the small de- 
gree nodes (compare the shaded regions on the left and 
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FIG. 2: (a) The mean queue size (S) for the small degree 
nodes (dashed line) and the hubs (solid line), (b) the mean 
queue size for the entire network (S) plotted as a function of 
a. Note that larger a correspond to larger nodal capacity. 

the right). Hence, we have far less control over the EE 
taking place on the small degree nodes. This is not the 
case for the hubs. The total number of EE on the entire 
network will depend on the events taking place on both 
the hubs and the small degree nodes. Thus, by adjust- 
ing the capacity or tuning a, we can manipulate the EE 
taking place on scale-free network. 

B. Excess load as queue size 

In any real network that encounters congestion-like sit- 
uation, the traffic in excess of the nodal capacity leads to 
a build-up of queue. For instance, clustering of vehicles 
at the traffic junctions or pending http requests to a web 
server or phone connections waiting to be serviced by cel- 
lular hubs are all examples of such a build-up. If a < 1/2, 
we should expect to see such queue and the number of 
walkers waiting in that hypothetical queue or buffer [] 
is an indicator of the severity of the extreme event. In 
models in which queue does not get cleared, under cer- 
tain conditions, jamming or a cessation of dynamics can 
take place. Phase transition to such a jamming state has 
been well studied ^14]. If Wi{t) represents the number 
of walkers on j-th node at times t = 1,2,...T, then the 
queue length is Qi{t) — {wi{t) ~~qi)9 (wiit) — qi). In this, 
unit step function 9{.) is used to ensure that Qi{t) = 
whenever Wi{t) < qi. We define the mean queue size (S) 
for a network with N nodes as 

NT N 
i=l t=l 1=1 

In this, (5) measures the mean number of walkers, in 
excess of nodal capacity, present in the queue per node 
at every time instant. Larger values of (S) imply more 
congestion-like scenario in the network. Now, we can ob- 
tain a broader picture of control over the extreme events 
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FIG. 3: Scaled mean number of extreme events {E/s,) / {Eo) 
with m = 4 and a = 1/2 as a function of strength of noise 
A. Random numbers were drawn from uniform and Gaussian 
distribution. 



as a function of a. Fig. [Hb) shows the mean queue size 
(5) obtained from random walk simulations with m — 4 
as a function of a. Evidently, (S) decreases with the in- 
crease in network capacity. Note that as argued before, 
a < 0.5 corresponds to a congestion-like situation and 
this is evident from the large values for (S). In contrast, 
for a > 0.5, (S) w indicating almost no extreme events 
on the network. Within the frame work of extreme events 
arising due to inherent fluctuations in the flux, it is im- 
possible to eliminate them altogether and hence (S) will 
never be exactly zero. However, the probability of ex- 
treme event occurrence and the value of (S) can be made 
arbitrarily small by tuning a. 

Even as we manipulate the overall congestion-like sce- 
narios on a network we explore the finer details. As Fig. 
[1] reveals, the extreme event probability on the hubs and 
small degree nodes have different dependence on degree 
k. This leads to an unequal apportioning of the burden 
of extreme events on hubs and small degree nodes of a 
network. We define a mean degree k such that we denote 
nodes with k > k as hubs and fc < A; as small degree 
nodes. Fig. [2Ia) show (S) for hubs (solid line) and small 
degree nodes (dashed line) separately. For a < 0.5, hubs 
display larger queue sizes as compared with the small 
degree nodes. In a reversal of roles, for a > 0.5, hubs 
display small queue sizes when compared with the small 
degree nodes. Quite interestingly, at a « 0.47, both the 
hubs and the small degree nodes equally share the bur- 
den of extreme events and queue sizes. This is indicated 
by the crossing of both the curves in Fig. [HJb). This is 
exactly the value of a at which the probability for ex- 
treme events becomes independent of degree k in Fig. 
[1] Irrespective of how the hubs and small degree nodes 
are defined, a « 0.47 would always remain the crossover 
point. Thus, except for a w 0.47, the network does not 
share the burden of extreme events in an egalitarian man- 
ner among the small degree nodes and the hubs. The 
practical implication is that the designing nodal capac- 
ity implicitly affects the 'spatial' distribution of extreme 
events over the nodes of a scale-free network. 
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III. EFFECT OF VARIABILITY IN NODAL 
CAPACITY 



In this section, we will take the capacity of i-th node 
to be Ci = Qi. The results in Sec. II implicitly assume 
that the capacity of i-th node d is related to its degree 
ki through Eq. [5] and that all the nodes with identi- 
cal degree have identical capacity. However, this con- 
dition is almost never satisfied in most of the real net- 
works [111, • In reality, the nodal capacity of i-th node 
does not follow any prescribed formula and indeed could 
be treated as a random variable drawn from a suitable 
probability distribution. Hence, we have for the time- 
independent capacity of i-th node 

a{A) = (/,) + {ma^ ± Aim) = C^{Q) ± A^m (7) 

with A being the strength of noise and a random num- 
ber. Evidently, if A = 0, this reduces to Eq. [SJ the 
threshold for extreme events. We compute the mean 
number of extreme events over the entire network, i.e, 

(Ea) = (l/A^T)E,r=iE.=i^KW -C^.(A)), scaled by 
{Eq). Thus, after simple manipulations, we get. 



4000r 



{Eo} J:^F^iC^m 



(8) 



By construction, this quantity is unity at A = 0. As sim- 
ulations in Fig. [3]reveal, {Ea)/{Eq) for m = 4 increases 
nonlinear ly as the strength of noise A increases. It agrees 
with the analytical result obtained using Eqs. [8] and ID 
Notice that Eq. [7] can be written as Ci(A) — {fi)-\- fhai 
where rh — m i A^i, i.e, every node has a different 
value of capacity indexed by m. We use the fact that 
the extreme event probabilities Fi scale with magnitude 
to obtain a leading order estimate of {Ea)/{Eq) 
l-|-ai A-|-a2A^, where ai and 02 depends on the degree of 
the node. This conclusion is independent of whether the 
random numbers were drawn from uniform or Gaus- 
sian distribution. 

Physically, this can be understood as follows. Suppose 
Ci (0) represents 4cr capacity for all the nodes of the net- 
work. For noise strength A > 0, some nodes will have 
capacity C_ such that C_ < Ci (0) while others will have 
a capacity C+ such that C+ > Ci(0). The nodes with 
C_ will encounter more extreme events and hence higher 
probability for the occurrence of extreme events (in com- 
parison with Ci(0) case). However for nodes with capac- 
ity C_|_ the extreme event occurrence probability is much 
smaller and leads to fewer events. Hence we notice a net 
increase in the number of extreme events on i-th node. 
As A increases, this effect translates into an increase in 
the mean number of extreme events over the entire net- 
work. In most of the real-life networks, the nodal capac- 
ity is typically proportional to the flux passing through 
the node 15| but is unlikely to adhere strictly to spec- 
ifications such as C(0) assumed in Eq. [71 Hence it is 
an appropriate guess that nodal capacity being propor- 
tional to its importance as measured by degree could be 
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FIG. 4: The ratio of change in mean number of extreme events 
S{E) on the network for SC change in capacity of the network 
plotted as a function of a. This is obtained from random walk 
simulations (solid circles) with m — 4. The inset shows the 
same data in log-log plot. 



its mean behaviour [111 but the actual capacities would 
be a random variable [16| clustered around this mean be- 
haviour. In all such cases, the random walk based model 
studied here predicts that the number of extreme events 
will increase with the variability in the nodal capacity. 
For other values of a too qualitatively similar results as 
in Fig. [3] are obtained. 



IV. CAPACITY ADDITION AND EXTREME 
EVENTS 

Given that the variability in the nodal capacity leads 
to an overall increase in the number of extreme events 
(as shown in Fig. [3]), we ask how we can reduce the 
mean number of extreme events on the network. This 
would be a desirable req uirement for a smooth function- 
ing of the networks [17]. As argued before, intuitively 
we expect mean number of extreme events to decrease 
if nodal capacity increases. We can increase capacity by 
increasing a. Further, to simplify the scenario, we take 
A = for this part. Figure S] shows simulation results 
for the change in the mean number of extreme events 
when capacity changes by one unit, i.e, 5{E)/5C as a 
function of a. For a < 0.5, the 6{E)/5C curve decays 
quickly. In this regime, capacity addition is accompanied 
by a large reduction in the number of extreme events. 
Hence, any capacity addition in this regime yields rich 
dividends in terms of cost-benefit ratio. However, be- 
yond a « 0.5, 6{E)/6C is vanishingly small indicating 
no significant change in the number of extreme events 
even when the capacities are increased. When a > 0.7, 
S{E) < 1, i.e, less than 1 extreme event on an aver- 
age. In this regime, capacity increase is not beneficial 
in alleviating the extreme events. Since capacity addi- 
tion (like building bridges, for instance) is generally an 
expensive proposition, we emphasize the following im- 
portant implication of this result. If the aim is to control 
the mean number of extreme events on a scale-free net- 
work even while maintaining a certain reasonable level 
of cost-benefit ratio, it is important to apriori estimate 
if we operate in the regime in which cost-benefit ratio is 
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favourable. The results obtained in this work provides 
insights into nodal capacity addition vis-a-vis reduction 
in EE. We emphasise that qualitatively similar results 
are obtained even when A > 0. 



V. CONCLUSIONS 

In this work, we have studied the extreme events in- 
duced by inherent fluctuations in the flux and taking 
place on scale-free networks. Given that such inherent 
fluctuations are unavoidable in any finite network, we 
show that it is possible to tune, either increase or de- 
crease, the probabilities for the occurrence of extremes 
in flux on a scale-free network using nodal capacity as 
a tunable parameter. This requires that the degree dis- 
tribution have a large support of at least two orders of 
magnitude (see x-axis in Fig. [Ij. This is not true of 
Erdos-Renyi networks with sharply peaked degree distri- 
butions and hence these results will not hold good for 
them. As we tune the nodal capacity, we show how var- 
ious parts of the network share the burden of extreme 
events and we discuss its implications. Further we study 
the effect of variability in the nodal capacity. We show 
that larger variability in the nodal capacity of scale-free 
networks leads to more extreme events. 



On the other hand, as we intuitively expect, increasing 
nodal capacity and hence the capacity of the entire net- 
work leads to decrease in the mean number of extreme 
events in the network. One significant result of this work 
is displayed in Fig. |4l This shows that increasing capac- 
ity beyond certain level does not lead to proportionate 
decrease in the incidences of extreme events. This has 
important implications for network design efforts. Since 
increasing the capacity of network is an expensive propo- 
sition in almost all the real-life situations, this work pro- 
vides an important theoretical benchmark to understand 
the limitations of capacity addition as a route to alleviate 
extreme events. In general, real life situations involving 
transport on scale- free networks are generally more com- 
plex than the random walk scenarios studied here 
However, this work sets the benchmark against which to 
understand extreme events arising from other types of 
transport processes on networks. Further, based on ev- 
idence from earlier studies 2], it is possible that even 
if random walk is replaced by a more intelligent routing 
algorithm, the results of this work might qualitatively 
remain valid. 

All the simulations were done on the computer clus- 
ter at USER Pune. ARS thanks DST for the financial 
support during the time this work was done. 
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